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We develop a thermodynamic description of particles held at a fixed surface potential. This 
system is of particular interest in view of the continuing controversy over the possibility of a fluid- 
fluid phase separation in aqueous colloidal suspensions with monovalent counterions. The condition 
of fixed surface potential allows in a natural way to account for the colloidal charge renormalization. 
In a first approach, we assess the importance of the so called "volume terms", and find that in 
the absence of salt, charge renormalization is sufficient to stabilize suspension against a fiuid-fiuid 
phase separation. Presence of salt, on the other hand, is found to lead to an instability. A very 
strong dependence on the approximations used, however, puts the reality of this phase transition 
in a serious doubt. To further understand the nature of the instability we next study a Jellium- 
like approximation, which does not lead to a phase separation and produces a relatively accurate 
analytical equation of state for a deionized suspensions of highly charged colloidal spheres. A critical 
analysis of various theories of strongly asymmetric electrolytes is presented to asses their reliability 
as compared to the Monte Carlo simulations. 

I. INTRODUCTION 

There is a long standing debate in the field of colloidal science concerning a possibility of phase separation in 
aqueous colloidal suspensions containing monovalent counterions P,|3,|3)0)ISIEI3>I1>0) ^> Qi ^1 ■ Charged 
colloids provide a particular challenge to the theorists. The long ranged Coulomb interaction and the extremely large 
asymmetry between the polyions and the counterions makes it very difficult to apply to this system the traditional 
methods of the liquid state theory. In particular, such a well established tool of condensed matter theorists as integral 
equations is found to be useless when applied to strongly asymmetric electrolytes at large couplings (low temperatures). 
For a wide range of parameters of physical interest the integral equations fail to even converge. The original hope of 
associating the lack of convergence with an underlying physical instability has proven to be unfounded |l5lJ l6 | . 

Recently, a linearized density functional theory has been used to study charged colloidal suspensions [ll 0| . The 
theory predicted a thermodynamic instability which manifested itself as a fluid-fluid phase separation. However, the 
underlying approximation of the theory put its conclusions in doubt 0, 0, • Specifically the linearization of the 
density functional lowers the statistical weight of the configurations in which the counterions are in a close vicinity 
of t he p olyions. This effect can be partially accounted for through the renormalization of the bare colloidal charge 
[TtLIISII . Unlike the bare charge, which can be very large, the effective (renormalized) charge is found to be bounded 
by the saturation value controlled — in a given solvent — by the colloidal size, temperature, and salt concentration. 
The renormalization of colloidal charge was argued to wash out the phase transition predicted by the linear theories 
0. Furthermore, numerical solutions of the full non-linear Poisson-Boltzmann (PB) equation inside a Wigner Seitz 
(WS) cell shows absence of any instability j^. lol IT^. On the other hand, linearization of the PB equation leads to 
a non-convex pressure, as a function of colloidal concentration, similar to the one observed in other linear theories. 
All these suggest that the phase instability predicted by the linear theories might be an artifact of the underlying 
approximations. 

To further explore these interesting points, we have investigated the thermodynamics of particles fixed at constant 
surface potential , 19.. ■20., .22.]. Relevance of the constant potential ensemble follows from the observation that unlike 
the bare charge, the effective charge of colloidal particles does not grow indefinitely but instead saturates. The 
saturation value is such that the effective electrostatic potential (ps in the vicinity of colloidal surface is Pqips ~ 4, 
where (3 = l/ksT and q is the elementary charge. Thus, although the colloidal charge can be very large, the potential 
near the colloidal surface, i.e. within the Debye length, does not increase beyond Pq(l>s ~ 4. This simple observation 
is sufficient to construct a consistent thermodynamic description of colloids with a state dependent effective charge. 
Since the WS cell description of colloidal suspension!^ does not lead to a fluid-fluid phase separation (see for a 
general argument), to further understand the mechanism of the instability observed within the linearized theories, we 
focus on this alternative treatment. 

We shall first (section In)) , study the effect of charge renormalization on the polyion-microion interaction free 
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energy, i.e. the volume term that appears when the original mixture of colloids, colons and counterions is mapped 
onto an effective one component system of dressed colloids P, • In section 11111 an alternative derivation of 
the thermodynamic equation of state based only on the far field considerations is presented. The corresponding 
pressure-density isotherms, in this case, do not exhibit criticality at any salt concentration. In order to assess the 
reliability of various approaches, in section Hvl wc compare the corresponding pressures to the results of the Monte 
Carlo simulations of Linse Q. We also consider the, recently proposed, symmetric Poisson-Boltzmann and the 
"boot-strap" Poisson-Boltzmann j^J theories. Conclusions are drawn in section Ivl 



II. ROLE OF VOLUME TERMS 

A. State dependent effective charges 

Consider a colloidal suspension at concentration pp, containing spherical polyions of charge —Zq and radius a 
in contact with a monovalent salt reservoir at concentration Cs- Now suppose one colloid is fixed at r = 0. In a 
continuum approximation consisting of smearing out the charge of other colloids and linearizing the Poisson-Boltzmann 
eauation[2^. the electrostatic potential at distance r from the center of colloid is 

0(r) = , B{x) = , (1) 

er (1 + X) 

where the inverse Debye length is 

K ^ V47rAB(p+ , (2) 

and the Bjerrum length is 

As = ^ . (3) 

The mean densities of colons and counterions inside the suspension are respectively p- and /9+. When pp — s- 0, 
^ = SttAs Cs. 

For highly charged polyions, Eq. ^ strongly overestimates the real electrostatic potential. However, it can be 
made consistent with the full non-linear PB, if instead of the bare charge Z an effective, renormalized, charge ^e// is 
used. The observation that for large surface potentials the electrostatics away from the colloidal surface is completely 
insensitive to the surface charge density allows for the "far-field" definition of the effective charge. Specifically, viewed 
from a distance larger than the Debye length, and provided that Ka > 1, the surface potential of a strongly charged 
colloidal particle appears to be Pq4>{a) ~ —4 0,|23,|22- Combining this with Eq. Q leads directly to [ial2^ 

Zeff = ^{1 + Ka). (4) 

We should note that this is the saturated value of the effective charge relevant for the highly charged colloidal particles. 
For weekly charged particles there is little or no charge renormalization. In the infinite dilution limit |23| of one colloid 
immersed in a 1:1 electrolyte of concentration Cg, the exact result for the saturation limit of ^e// is 

^e// = ^ (2 ■ 

This expression -valid up to corrections of order (ko)"^, which turn out to be quite small as soon as ko > 1, is very 
close to the approximation |0J, derived from the matching procedure detailed in [2Q; 22]. 

The effective charge depends strongly on the electrolyte concentration inside the suspension. Salt screens the 
electrostatic interactions between the counterions and the polyions and leads to an increase in the colloidal effective 
charge. While the effective charge of colloidal particles is a strongly state dependent function, the effective surface 
potential is not. Therefore, inside suspension, colloids behave as if their surface potential was effectively fixed. It is 
interesting, therefore, to study the thermodynamics of colloidal particles at fixed surface potential. 
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B. Thermodynamics of particles at fixed surface potential 

The change of the thermodynamic ensemble from the constant charge to the constant surface potential allows us in 
a natural way to explore the role of charge renormalization in the framework of a linear theory. The surface potential 
is related to the effective colloidal charge through the Eq. 

</5 = -M{a) = /-f^^ N ■ (6) 
a(l + Ka) 

As mentioned above, the reduced surface potential within the non-linear Poisson-Boltzmann theory is found to 
saturate at (p = 4. For the sake of generality we shall, however, keep its value arbitrary. In the subsequent analysis, 
Z will refer to the saturation value of the colloidal effective charge. An implicit assumption is therefore that the bare 
charge largely exceeds the effective one. 

In the simplest approximation, the Helmholtz free energy of the suspension is a sum of entropic and electrostatic 
contributions 

PF = N+[liip+A^ - 1] + Ar_[inp_A3 - 1] + /JF'^' , (7) 

where A is the de Broglie thermal wavelength, and N-^- and N- refer to the number of counterions and colons inside 
the suspension. 

The electrostatic free energy results from the polyion-microion, microion-microion and the polyion-polyion in- 
teractions. For suspensions containing monovalent counterions the polyion-counterion interaction is the dominant 
contribution and will be the only one kept in the present exposition (we come back to this point in section llV)l . We 
find[23 

-2a{l + Ka) ■ 

This expression can be obtained through the usual Debye charging process in which all the particles are simultaneously 
charged from to the their final chargej^^. Alternatively a surface charging process, at constant Debye length can 
be employed|23l- It is noteworthy that in the salt free case, we recover precisely the volume term obtained in Refs 

m 

If the suspension is in contact with a salt reservoir of chemical potential /is , the effective charge of colloidal particles, 
as well as the number of counterions and colons, is determined by the minimum of the grand potential function 



ip{l + Ka) 



7(7V+ - iV_ - ZNp) . (9) 



In this equation, uj and 7 are the Lagrange multipliers: 7 ensures the charge neutrality of the system, while uj 
enforces the relationship between the surface potential and the effective charge, Eq. 0. In the biophysics literature 
7 is known as the Donnan potential. It results from the inability of macroions to diffuse through a semi-permeable 
membrane. This is precisely the situation that we have in mind, while the microions are assumed to be in a free 
exchange with the reservoir, the polyions are confined to the interior of suspension. This restriction on the polyion 
mobility results in a potential difference 7, between the bulk of suspension and the reservoir. The colloids are then 
held at potential ip + "f with respect to the reservoir, or equivalently at potential if with respect to the bulk of the 
suspension. Within the WS cell model, a similar prescription of constraining the potential difference between the 
colloidal surface and the outer (reservoir) boundary has been shown to yield a surprisingly good agreement with the 
full non- linear Poisson-Boltzmann eauation|2Clj|. 

Minimizing the grand potential with respect to 7V+, N- and Z, we find 



dpn _ d(3F 



2'KLp\Bapp 



T7 = 0, (10) 



dl3n df3F 



dZ 



dZ 



"As" 

a 



ipNp = 



(11) 



and the charge neutrality condition reads 



ipa{l + Ka) 



(12) 
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Noting that 



ipNp , (13) 



dZ 

Eq. simplifies to 

cj = . (14) 

Eliminating the Lagrange multiplier between Eqs. H1U|) . we are left with two equations which govern the concentrations 
of counterions and coions inside suspension, 

P+P- = exp I ^ 1 (15) 

and 

Lpa{l + Ko) 

P+- P- ^ T Pp ■ (16) 



C. The equation of state 



The osmotic pressure inside suspension of colloids at fixed surface potential is 



Np,fj.s,(jj,f,ip 



(17) 



It is important to keep in mind that as Q changes with volume, the number of coions, counterions, as well as the 
charge of colloidal particles are all varying. This is the reason for writing the total derivative in the expression H17[l 



dV 



dn dN+ dn dN_ 



dN+ dV dN_ dV 



dn dz 

dZdV 



dn 

dV' 



Recalling that at the thermodynamic equilibrium 

dn 



dn 



dN+ dN. 



dn 

dZ 



expression for pressure simplifies to 



P 



dn 

dV 



Np,N+,N-,Z,fj,s,uj,^,ip 



dF 
"dV 



N„,N+,N^.Z 



(18) 



(19) 



(20) 



Equation (|2()|l beautifully illustrates the thermodynamic principle of ensemble equivalence. The functional form of the 
pressure is the same weather the calculation is done in the fixed potential ensemble using the grand potential function 
n, or in the fixed colloidal charge ensemble using the Helmholtz free energy F. We stress that simply inserting 
Z(,ff{V) with its state dependence into F and then differentiating it with respect to volume will lead to an incorrect 
result. If the Helmholtz free energy is used, the variation must be performed at fixed colloidal charge. 

In general it can be very difficult to find a suitable thermodynamic potential for a constrained system. The 
calculation of pressure, on the other hand, can be done very straightforwardly using the constant Z ensemble, and 
enforcing the constraint a posteriori. Evaluating the partial derivative in Eq. I|20(l . the osmotic pressure inside the 
suspension takes a particularly simple form 



f3P = P+ 



1 a 



■—f napp , 



(21) 



where the concentrations of coions and counterions are determined from Eqs. H15|l and ()16|1 . 
In the special case of vanishing salt concentration, Eqs. (|15|l and Hlt)|) simplify to 



(22) 
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and the ratio of colloidal size to Debye length is 



^+1^3w(4 + 3w) 



where 77 = AirppO? /'i is the macroion volume fraction. For salt-free suspensions pressure becomes 

2 f 1 
47rABa [iP — irjip 1 + na (pKa 



(23) 



(24) 



(25) 

li if < A the osmotic pressure is a convex up function of colloidal density. For real colloids with = 4, the pressure 
is a linear function of colloidal density, PP = Appa/Xs- For surface potentials strictly above 4, a thermodynamic 
instability appears. It is very curious to note that the instability sets in precisely at = 4, which is the saturation 
value for the surface potential obtained within the non-linear Poisson-Boltzmann theory. Renormalization of the 
electrostatic free energy is, therefore, sufficient to stabilize a salt free real colloidal suspension [ip — 4) against a 
fluid- fluid phase separation '7'|. However, the fact that the critical surface potential is precisely equal to the saturation 
value of the non-linear theory, suggests that the approach, most likely, is very sensitive to the approximations made. 
Furthermore, a relatively small amount of salt destabilizes suspension even when (/? = 4, see Fig. ^ The critical 
salt concentration for highly charged colloids of radius a = 1000 A, is c* ~ 10~* M, which corresponds to K*a ~ 3.3. 
Since there is no explicit polyion-polyion nor niicroion-microion interaction, the instability is completely driven by the 
polyion-counterion correlations. In Fig. |2]we show the effective charge Z resulting from our approach, as a function 
of volume fraction. A good agreement with the Poisson-Boltzmann cell model is found. 
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FIG. 1: Pressure isotherms for different reservoir salt concentrations (k^ — SttAsCs, Lp — A) 

It is important to stress that even a minor modification of the approximations employed may have a dramatic effect 
on the predicted phase instability. One may wish for instance to use Eq. ij^l for the effective charge instead of Eq. 
Q. It is important to note, however, that such a modification means that colloid is no longer held at fixed potential. 
Therefore, the grand potential SI, as written in 0, can no longer be used. However, we can compute the functional 
dependence of P on Z , pp and Cg by differentiating the Helmholtz free energy F with respect to volume at constant 
Z, and enforcing the constraint Z = aiAna -I- 6)/ As a posteriori. Following this route, we recover the same equation 
of state as before [i.e. Eq. (|25|l in the salt free case but with now a different salt dependence of Z] and the critical 
salt concentration above which the instability sets in (see Figure 13, decreases by a factor of four to K*a ~ 1.77. 



III. AN ALTERNATIVE APPROACH : A JELLIUM APPROXIMATION 



In the previous section, we have found a fairly accurate expression for the effective colloidal charge at saturation 
and used it to renormalize the electrostatic free energy. It is important to remember, however, that the effective 
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FIG. 2; Comparison between the effective charge (sohd hne) found using our variational approach, to that obtained within the 
Poisson-Bohzmann cell theory (circles), following the prescription proposed by Alexander et al. l34j 
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FIG. 3: Pressure isotherms for different reservoir salt concentrations, making use of Eq. Q instead of Eq. The situation 
is now different from that of the previous constant <f ensemble. Working is the constant Z ensemble with the Helmholtz free 
energy (|7|l nevertheless allows to compute the pressure. 



charge, is by definition related to the "far field" asymptotic properties of the electrostatic potential. Its use for the 
renormalization of electrostatic free energy is, therefore, questionable, since not only far field but also the near field 
properties of the electrostatic potential may be relevant. We now turn our attention to a simple approach which relies 
only on the far field features of the electrostatic potential to obtain the equation of state. As within the WS cell 
picture, use will be made of contact theorem, which relates the osmotic pressure to the concentration of counterions 
in the region where the electric field is zero|2ll|. 

We now reconsider the approach put forward at the beginning of section^ Consider one colloidal particle fixed at 
the origin of coordinate system. As before the charge of other microions is uniformly smeared throughout the solution. 
On the other hand, positions of counterions are strongly correlated with those of colloids. The system then forms 
a jellium, where the electrostatic potential far from the colloid 0oo (bulk) differs from that in the reservoir (chosen 
to vanish). The solution of the linearized Poisson-Boltzmann equation for Scf) — (fi — (j)oo is again given by Eq. l|T]l. 
where the screening length is related to the bulk salt concentration, = 47rAB(/5+(oo) + p_(oo)). In the spirit of the 
previous discussion, the colloidal particles are held at constant surface potential = 4 with respect to the bulk (poo, 
which again imposes Z — 4a(l + Ka)/XB- The counterions and colons are distributed inside the jellium in accordance 
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As/a 


pP/n, MC 


PP/n, Z = ^eff 


PP/n, Z = Zbarc 


Zeal Zhmc 


0.022 


0.98 


0.99 


0.99 


1.0 


0.044 


0.95 


0.96 


0.97 


0.99 


0.089 


0.89 


0.88 


0.92 


0.95 


0.178 


0.71 


0.68 


0.78 


0.82 


0.356 


0.45 


0.41 


0.46 


0.53 


0.712 


0.26 


0.21 


-0.32 


0.29 



TABLE L Equation of state 12511 with (third column) and without (fourth column) renormalization of volume terms, as a 
function of electrostatic coupling As/a, for a packing fraction 77 = 0.00125. The quantity n denotes the mean total density 
of counterions n = /9p.^bare- The MC data (second column) are taken from reference 6]. Since the previous parameters do 
not correspond to the saturation regime of effective charges but only approach it, we have used the effective charge given by 
Alexander's prescription [l^ to compute the pressure from Eq. 1)25^ in the third column. The corresponding ratio .^/^barc is 
indicated in the last column 



with the Boltzmann distribution 



P+{r) = Cs 



-I3q<t>{r) 



(26) 



(27) 



where (j>{r) is the local electrostatic potential with respect to the reservoir. Taking the product of Eas. (|26(l and 
(|27|l we find the famihar condition for Donnan equilibrium /9+(r)p_(r) = c^. The electro neutrality constraint 
p+(oo) = p-{oo) + Zpp, closes the problem j25l|. The total concentration of niicroions inside the suspension is then 
related to their concentration inside the salt reservoir through 



[naf = {Ksaf + [12ry(l + Ka)f 



(28) 



The osmotic pressure, within the non-linear PB theory, is determined from the concentration of microions in the 
region where the electric field is zero. Since the electrostatic potential decays exponentially with r, the electric field 
vanishes when r — > co. Within the jellium approximation the osmotic pressure then takes a particularly simple form, 



One can show that solution of Eq. (|28|l obeys the inequality 



47rAf 



Ok 
dr] 



> 



(29) 



(30) 



which ensures that the compressibility is always positive and that suspension is stable against the phase separation. 
In the absence of salt, we obtain a simple analytic expression 



477X30^ f3P = 1277 1 + 677-1- Vl2?7(l + 377) 



(31) 



that will be tested against experimental data in section iTVl 



IV. DISCUSSION 



A. Comparison with the Monte Carlo simulations 



The above analysis shows that different routes to thermodynamic pressure lead to very different results. In order to 
decide which route is the most reliable, a comparison with "exact results" is welcome. As a benchmark, we can use the 
Monte Carlo (MC) pressure data of Linse for salt-free asymmetric electrolytes consisting of highly charged spherical 
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FIG. 4: Monte Carlo pressures (dots) compared to those obtained within PB cell model (continuous curve), boot strap PB [Tl| 
(crosses) and symmetric PB theories [lol (triangles, corresponding to virial, charging and compressibility routes). The packing 
fraction is 77 = 0.00125 and n — Pp^barc- in table|l] the charge asymmetry polyion/counterion is .^barc — 40. 



macroions and point counterions. At high electrostatic couplings, this system exhibits an instability and separates 
into two coexisting phases of different electrolyte concentration. We shall argue that a minimum requirement for a 
reliable theory of phase behavior is its ability to reproduce reasonably accurately the MC equation, at least up to the 
transition point. This appears to be a stringent test and a necessary condition to trust any instability that a theory 
might predict. 

We first test in table^the equation of state (|25|l for a charge asymmetry Zbaro = 40 between colloids and counterions. 
It is evident that renormalization of colloidal charge significantly improves upon complete neglect of non linearities. 
The latter approach consists in considering Z ~ ^barc and severely fails at high Xb/o,, giving negative pressures. 
Our renormalized volume term captures the main effect of nonlinearities, but the agreement with MC, even if decent 
in view of the simplicity of the approach is nevertheless only qualitative, and does not reach the level of accuracy 
required to discuss phase stability. 

The jellium equation of state derived in section ITTll onlv holds for saturated effective charges, i.e. in a regime of 
coupling that the Monte Carlo simulations, so far, have not reached (which corresponds to a very high bare charge 
with a large separation of scales between Bjcrrum length and colloid radius, see below). We therefore directly turn to 
the comparison of the relative performances of the Poisson-Boltzmann cell model, symmetric PB 1Q| and boot-strap 
PB (see for details), with respect to Monte Carlo data (see Figures 0] and [3) . 

To produce these figures (providing a similar comparison as Table , we have chosen the lowest and the highest 
packing fractions investigated by Linse in lol. The striking feature revealed by Figs 0] and is the remarkably good 
agreement between the PB cell pressures l2a. l27l and the MC simulations, even at = 0.00125 where the cell model 
could have been anticipated to fail (see also [23). The only competitive approach at 77 = 0.08 seems to be the boot 
strap PB theory but this theory severely fails for low volume fractions, see Fig0| necessary to study colloidal 
phase-stability. At this volume fraction, the simple treatment of section ^ [Eq. (|25|l ] , provides a better equation of 
state than the boot strap PB or symmetric PB (see Tableland Fig. QJ- Within the PB cell, polyion-polyion as well 
as counterion-counterion correlations are discarded ; Figs 0] and |5l show that as far as the pressure is concerned, these 
contributions are small or negligible, for the parameters investigated, even at the highest couplings. This justifies 
their neglect in our analytical treatment. 

Given the accuracy of the PB cell model, we are now in the position to assess the quality of the jellium approximation 
of section Unl The corresponding pressures are compared in Fig. |Hlwith their PB cell counterparts for highly charged 
colloids (where the effective charge saturates to its upper threshold), both with and without added salt. The simple 
analytical expression for salt free suspensions is found to be in good agreement with the PB data. Unfortunately 
the agreement deteriorates when KgO > 1 (see the inset). 
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FIG. 6: Comparison between PB cell pressures and those for the jellium model, without added salt [in the latter situation, the 
equation of state is given by expression Inset: same when the suspension is dialyzed against a salt reservoir such that 

Ksffl = 2.6 (the quantity P considered is the osmotic pressure, i.e. the reservoir contribution has been subtracted). 



B. Relevance for colloidal suspensions 



At this point, we must conclude that the PB theory, even restricted to the cell, is superior to the competing 
approaches for aqueous suspensions with monovalent counterions. At high electrostatic couplings, corresponding to 
multivalent counterions in water, the MC simulations of Linse 6] find an instability. This transition, however, has 
nothing to do with the volume terms, but is the result of strong correlations between the double layers of colloidal 
particlesji^ which produce attraction between like-charged colloids at sufficiently short separations[23 133- So far, 
this attraction has not been properly included in any of the thermodynamic theories of colloidal stability. 

Validity of PB theory. To quantify the range of validity of the PB theory one may construct a dimensionless 
parameter F^c characterizing the importance of microions correlations, discarded within the PB theory. For monovalent 
microions. Fee oc /3g^/(€e) where £ is the characteristic mean distance between the microions in the double layer. If 
the number of condensed counterions is such as to almost completely neutralize the colloidal charge, which is the case 
for strongly charged colloids, i ~ a/y/Zhaic, and Fee becomes [23 . l29l| 



(32) 
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FIG. 7: Same as figure El including a comparison with osmotic pressures of deionized bromopolystyrene particles (shown by 
diamonds, from ref. ^33]). The curve labeled "jellium" corresponds to the analytical prediction l|31|l . No adjustable parameters 
have been used. 

When Tec exceeds unity, PB theory is expected to break down. The value Tec — 2 has been reported to correspond 
to the instability threshold |^ , which and has been observed in the simulations of Linse . The field theoretic 
treatment of Netz also corroborates this conclusion For particles with Zbaro = 40, Fee — 2 corresponds to 

Xs/a — 1.1. Thus the PB theory can be expected to work quite well up to very high surface charge concentrations. 
Indeed, comparing the predictions of the PB cell model to the MC simulations, an excellent agreement is observed 
up to As/a ~ 0.7 where the MC data stop, rather close to the expected point of instability As/a ~ 1.1, see Figs. 01 
and [SI 

Validity of the saturation picture within the PB. The constant potential approach used in sections ^ and IIIII relies 
on the phenomenon of effective charge saturation exhibited by the PB theory, when Zbare is large enough. The 
saturation occurs when the electrostatic energy of the condensed counterion is significantly larger than /csT. This 
can be characterized by a dimensionless parameter 

-^bare -^barcAi^ /qqN 

t sat = = JZ ' T- (.O-jj 

When Fsat becomes larger than 1 |20ll23| . linearized theory fails and charge renormalization becomes important. 

We must stress that large values of Fsat are fully compatible with small values of of Fee- Specifically, for any 
a ^ As, there exists a range of bare colloidal charges Zbaro, such that Fee < 2 < Fgat- For these values of Zbarc, 
the polyion-microion interaction is sufficiently strong that linearized theories, without charge renormalization, will 
certainly fail. On the other hand the counterion-counterion correlations are sufficiently weak, so that the PB theory is 
still applicable. To illustrate this point, we compare in Fig. [Tithe pressure obtained at saturation within the PB cell 
model calculation (formally Zbare oo), to that measured experimentally by Reus et al, under the conditions very 
close to complete deionization (no salt). The agreement with PB theory had already been mentioned in (20ll33l |. but 
provides an illustration of the saturation phenomenon in real suspensions. It also shows that despite its simplicity, 
Eq. (|31l) is fairly accurate. 

It is important to stress, however, that for any colloidal size a, there is a maximum value of .^bare above which 
the PB theory fails. In practice, however, this break down of the PB equation never occurs for aqueous suspensions 
containing only monovalent counterions. The reason for this is that the hydrated ionic size provides a lower cutoff for 
the length £ which appears in Fee, i-e. i > d, where d « 4 A, is the hydrated ionic diameter. For aqueous suspension 
with monovalent counterions Fee is, therefore, always less than Xs/d, so that Fee < 2. For monovalent counterions in 
water, the PB theory, therefore, should apply without any restriction. Of course, this pleasant situation changes as 
soon as multivalent salt is added to suspension. In this case, hydrated ionic size is no longer sufficient to restrict the 
strength of microion-microion correlations and Fee > 2. Under these conditions the PB theory will no longer apply 
and a more sophisticated approach must be used. 
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V. SUMMARY AND CONCLUSIONS 

We have proposed (section |nj a linear theory to investigate the phase behavior of colloidal suspensions. The non- 
linear effects are partially taken into account through the postulate that highly charged polyions behave as if they 
were constant potential objects. The effective charge of colloidal particles is, therefore, a state dependent function. 
Our first goal was to develop a consistent thermodynamic approach for such state dependent charges. The results 
found in section^ rely on a simple form for the electrostatic free energy (volume term), resulting from the polyion- 
microion interactions calculated using the linearized PB equation, for spherical colloids. The approach could be easily 
generalized to the case of cylindrical macroions. The critical behavior predicted is, however, spurious, which may 
be attributed to the simplicity of the volume term used, and/or the difficulty of renormalizing such terms. The 
jcllium-like model of section ITTll provides a more reliable route, and allows to obtain analytically a simple equation of 
state for highly charged colloids in the salt free limit, see Fig. Unfortunately it is difficulty to see how this kind of 
approach can be extended to account for the polyion-polyion interactions. 

^From our analysis, we conclude that for Tec — [■^barc/(47r)]^/^AB/a < 2, Poisson-Boltzmann approach (PB), even 
restricted to the cell model, leads to more accurate predictions for the thermodynamic functions than the competing 
theories. This is quite remarkable, since it is by far the simplest (see the appendix of js^l for a "ready-to-use" 
implementation of PB cell model). However, there is clearly a need to go beyond the PB theory when dealing with 
the multivalent counterions, since it is the counterion-counterion correlations that drive a phase instability for Tec > 2. 
Inclusion of these effects in a theoretical approach is a difficult task, since they have little signature on the pressure 
data up to the electrostatic coupling where suddenly they destabilize the system. 

This work was supported in part by the Brazilian agencies CNPq and FAPERGS and by the french CNRS. Y. Levin 
acknowledges the hospitality of the Theoretical Physics Laboratory in Orsay where part of this work was performed. 
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